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TRANSIENT  PROCESSES  IN  AN  INDUCTOR  SYSTEM  CONSISTING  OF  A  FLAT  COIL 
AND  A  MULTILAYER  CONDUCTING  MEDIUM 

V.  N.  Bondaletov,  V.  P.  Gel'yetov,  Ye.  N.  Chernov 

Introduct ion 

Flat  inductors  placed  over  a  conduct ing  medium  are  used  in  a  whole 
series  of  electrical  engineering  devices  (during  electroinduct ive  flaw 
detection,  induction  heating,  magnetic  pulse  treatment  of  metals, 
etc .  ) . 

Steady  electromagnetic  processes  in  these  fixed  inductor  systems 
are  usually  calculated  in  the  approximat ion  of  a  plane  electromagnetic 
wave.  The  stationary  electromagnetic  field  1  of  an  annular  coil  with 
a  sinusoidal  current  placed  over  amultilayer  medium  was  also  found  by  solv¬ 
ing  the  heterogeneous  Helmholtz  equation  for  the  vector  potential  of 
the  magnetic  field. 

The  transient  electromagnetic  processes  which  take  place  during 
the  discharge  of  a  capacitve  store  to  an  inductor,  mainly  characteris¬ 
tic  of  magnetic  pulse  deformation  of  metals  and  dynamic  induction 
acceleration  of  conductors,  are  calculated  in  the  approximation  of  a 
plane  electromagnet ic  wave  2.  However,  the  assumption  of  the  presence 
of  a  single  field  strength  component  does  not  always  turn  out  to  be 
true  in  practice.  Therefore,  in  this  work  we  calculate  the  transient 
processes  for  the  more  general  case  -  a  flat  disk  coil  with  a  current 


placed  over  a  multilayer  conducting  medium.  We  will  use  a  coil  with 
a  surface  current  uniformly  distributed  over  the  dimensions  and  placed 
horizontally  over  an  unbounded  medium  (Fig.  1)  as  the  calculation 
model . 


The  field  of  a  real  inductor 
can  then  be  found  by  adding  the  fields 
of  the  indi.idual  disks. 

1.  Calculation  Method 

The  initial  equations  for  solv¬ 
ing  any  electrodynamic  problems  are 
the  Maxwell  equations,  which,  dis¬ 


regarding  the  bias  currents,  we  will 


write  as  follows  for  isotropic  media: 


*/em% 

rocc  H  , 


(  1) 


where  J  is  the  specific  conductivity  of  the  medium;  Jem  is  the 
density  of  the  extraneous  currents,  which  are  controlled  by  an  extern¬ 
al  source. 


We  will  make  the  calculation  by  the  successive  application  of 
the  integral  equations. 

We  will  use  the  Laplace  transformation  to  eliminate  differentia¬ 
tion  through  time  t. 

Then 


(2) 


\ 


where  A(p) 
Laplace);  n 


is  the  expression  for  the  vector  potential  (according  to 
is  the  relative  magnetic  permeability  of  the  medium. 


Solving  system  of  equations  (2)  analogously  to  [3],  we  arrive  at 
a  heterogeneous  Helmholtz  equation  for  expressing  the  vector  potential 


where 


In  the  future  we  will  bear  in  mind  that  all  of  the  transforma¬ 
tions  are  made  with  transforms  of  the  functions,  without  stipulating 
this  specially. 


We  will  use  a  cylindrical  coordinate  system  with  the  z-axis 
directed  normal  to  the  surfaces  of  the  layers  and  coinciding  with  the 
axis  of  the  disk  for  the  calculation;  then  (3)  changes  into  a  second- 
order  equation  in  partial  derivatives  relative  to  the  only  (because 
of  the  axial  symmetry)  component  /4yt 


We  will  eliminate  the  differential  operations  on  the  coordinate 
r  according  to  (3]  using  the  Hankel  integral  transformation  with  a 
kernel  in  the  form  of  the  first-order  Bessel  function 

Xs  JrZMAfop,  z)dr .  ( 5 ) 

.  0 

where  is  the  transformed  vector  potential;  r  is  the  weighting  func¬ 
tion;  ZN  is  the  kernel;  A  is  the  transformation  variable. 

Then  the  problem  in  question  is  reduced  to  an  ordinary  second- 
order  differential  equation  relative  to  the  unknown  twice-transformed 
f  unct ion 

(6) 


3 


where 


13 ' 


he  root  with  the  positive  real  part 


The  general  solution  for  nonhomogeneous  equation  (6)  is  known 


**>(-? 

**v(-p)(e*ffcrn  , 


where  £  is  the  integration  variable  along  direction  z;  B  and  C  are 
arbitrary  constants  determined  from  the  boundary  conditions. 

As  we  know,  the  conditions  of  continuity  of  the  vector  potential 
and  the  tangential  component  of  the  magnetic  field  strength  must  be 
satisfied  on  the  interface  of  the  media: 


32  *z 


with  Z  =  Zo. 


After  finding  the  transformed  function  ,  the  unknown 

vector  potential  can  be  found  by  using  the  successively  inverse 
transformations  of  Hankel: 


and  Laplace: 


y 


2.  Field  of  Disk  Coil  Located  Over  Two-layer  Medium,  Plate, 

Half -Space 

We  will  find  the  expression  for  the  vector  potential  in  the  suf¬ 
ficiently  general  case  when  a  disk  coil  is  located  over  a  layer  of 
thickness  d  with  conductivity  /j  lying  on  conducting  half-space  (JZ  J. 

We  will  set  pi  =  M2  =  P 3  =  1  >  which  is  valid  for  nonferromagnetic 


1 


media. 


For  the  upper  half-space  (z>0),  considering  that  .  •  0 
,  from  (7)  we  obtain 

X*  {jtfvMfaft*  *y>M)d$)+  e*fi(-\z)(cf  * /fa  •*>(*$)</ $)J \ 


In  view  of  the  absence  of  extraneous  currents,  we  will  have  the 
following  for  the  conducting  layer  (Z<0): 


4  'if/?  'rfc  *b 


For  the  lower  half-space 


4  • 


since  C3  =  0,  for  when  ,  the  field  should  be  bounded. 

The  field  should  also  be  bounded  as  ;  therefore,  from 

(11)  we  can  determine: 


4m/fa«p£Afy4. 


We  will  consider  the  extraneous  current  density  to  be  independent 
of  r  and  equal  to 


when  Z~6ir,*r>0 
/or,  faff  when  Z4h't  r>J}f 


Then,  with  consideration  of  (16),  we  obtain  the  following  for 
transformed  value  : 


fa- z  fcr*)  fa 
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where  is  the  Dirac  delta  function 

F(A)=  f rj/krjjr- 

Substituting  (16)  in  (14),  we  will  find  the  constant  Bi 

4  =f],/:M*xP ('* *  },£(*) exp  fat). 


(17) 


We  will  use  boundary  conditions  (8)  to  find  constants  B2,  B3, 

C-),  C2.  Substituting  the  values  of  the  vector  potentials  from  (11-13) 
and  their  derivatives,  on  the  interface  of  the  media  we  will  have: 


fat  e*flf  ?j  d)*cx e*P(?2 
4  up('h  rf'Pi  expffy  (/)*  4  3  d)- 


(18) 


Solving  this  system,  we  will  obtain: 


4- fa  .  s,-  f'l'Le.pffrfS,’/,*,  ; 

Substituting  these  values,  and  also  (17),  in  (11-13),  we  will 
obtain  the  expression  for  the  transformed  vector  potentials: 


(  19) 


exp(.?>l)crp( . 


(20) 


(21  ) 


(22) 
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Expressions  (20-22)  can  be  used  as  the  initial  expressions,  e.g., 
for  calculating  the  pressure  during  magnetic  pulse  deformation  of 
thin-walled  articles  in  a  metal  die. 


If  a  conducting  plate  of  thickness  d  is  located  on  an  insulating 
base,  then  fjm&  »  and  the  expressions  for  the  transformed 

vector  potentials  assume  the  form: 


(23) 


(24) 


X-u  7  r(\ } t -Kty 


(25) 


In  the  limiting  case  (  •  ),  which  corresponds  to  the  system 

of  a  disk  coil  over  a  half-space,  expressions  (23-25)  are  considerably 
simplified: 

/£-  e*p(-**  -*A},  (  26 ) 

j7T~e*P(g*-d)  •  (  27  ) 


3.  Transient  Processes  in  the  System  "Disk  Coil  -  Conducting 
Half-Space"  When  a  Sinusoidal  Decaying  Current  is  Connected 
to  the  Source 


The  connection  of  an  inductor  system  to  a  sinusoidal  decaying 
current  source  is  typical  of  induction  casting  of  massive  conductors, 
"  well  as  of  magnetic  pulse  treatment  of  metals.  In  this  case,  the 
current  density 


fi*}.  Ale, 


(28) 


i 
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where  a  is  the  attenuation  decrement  and  w  is  the  circular  current 
f  requency . 

Substituting  the  value  /*{/>}  in  formulae  (26),  (27)  and  using 
inverse  transformations  (9),  (10),  we  obtain  the  following  expressions 
for  the  original  vector  potentials: 

V/-  * 

.  .  -  v  (29) 

"  ftje pfa/>-  Ai)d*dp  --W  ^FL^c/pefr .. 

<30) 

W" 

^xp(-Ah^z)dpdX  - r) l */“* ■ 

We  will  introduce  the  variable  ;  then  for  /y^  ,  according  to 

[5,  P-  657],  we  will  have: 

Ffa)s  rfJr'J/hr'Jctr  (3D 

where  r'*fj  Wn)'JXn)W’H'tHZrn)  l 


H  is  the  Struve  function. 

Integrating  (29),  (3)  on  a  complex  plane,  after  simple  transform¬ 
ations,  with  considerat  ion  of  (3D*  we  obtain  the  following  expres¬ 
sions  for  the  vector  potentials: 

‘  TiX far’)  fyt)  wC-**  *d~  expf'*t)co$  r. 

/ft*  ('& gflr) dnj) 


8 


\yf* 

( 


where 


*'-i> 


1  s 


*-■ 


the  depth 

A 


of  penet ra t i on ; 


Having  the  expressions  for  the  vector  potentials,  we  can  determ¬ 
ine  all  the  values  which  characterize  the  electromagnetic  field  in 
the  system  in  question. 


4.  Magnetic  Field  Pressure  on  a  Massive  Conductor 


We  will  calculate  the  currents  in  a  massive  conductor  and  the 
magnetic  field  pressure.  This  problem  arises,  f-r  example,  when 
s t ra i gh t en i ng  large  articles  by  the  magnetic  pulse  method.  The  mag¬ 
netic  field  pressure  is  usually  determined  from  the  volume  Jensity  of 
the  pondermotive  forces: 


(  31*) 


where  J  is  the  volume  density  of  the  currents  in  the  conductor;  5  is 
the  magnetic  field  induction. 


The  pressure  along  the  z-axis  (Fig.  1)  is  obtained  as 

* 

%  &  , 


(  35) 


where 


Considering  that  during  magnetic  pulse  deformation  of  metals, 
■£*«  /  and  the  current  damping  is  small,  after  simplification,  we 

n 

can  obtain  the  following  from  (33): 

\ 

^  (37) 


**P/'-2*Qso,*r , 


(38) 


where 


Fig.  2.  Distribution 
of  pressure  on  surface 
( h  *  =  0.1). 


Thus,  the  dependence  of  the  pressure  on  time  is  the  same  as  in 
the  case  of  a  plane  electromagnetic  wave.  However,  the  pressure  dis¬ 
tribution  on  the  surface  does  not  remain  constant,  since 


(39) 


The  dependences  Bo(r)  and  Po(r),  obtained  by  numerical  integration  for 
the  case  h'  =  0.1,  are  given  in  Fig.  2. 


Conclusions 


The  solution  of  the  Maxwell  equation  for  the  "flat  disk  coil  - 
multilayer  conducting  medium"  inductor  system  with  an  arbitrary  depend¬ 
ence  of  the  exciting  current  on  time,  obtained  by  the  successive 
application  of  the  integral  Laplace  and  Hankel  transformations,  made 
it  possible  to  obtain  the  general  expressions  for  the  vector  potentials 
in  different  media.  The  expressions  for  the  vector  potentials  were 
also  obtained  with  the  connection  of  a  "coil  -  conducting  half-space" 
system  to  a  decaying  sinusoidal  current  source. 

The  inverse  transformat  ions  were  made  by  numerical  integration. 
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